Manipulating the rotational properties of a two-component Bose gas 
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A rotating, two-component Bose-Einstein condensate is shown to exhibit vortices of multiple 
quantization, which are possible due to the interatomic interactions between the two species. Also, 
persistent currents are absent in this system. Finally, the order parameter has a very simple structure 
for a range of angular momenta. 
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When a superfluid is set into rotation, it demonstrates 
many fascinating phenomena, such as quantized vortex 
states and persistent flow [l[ . The studies of rotational 
properties of superfluids originated some decades ago, 
mostly in connection with liquid Helium, nuclei, and neu- 
tron stars. More recently, similar properties have also 
been studied extensively in cold gases of trapped atoms. 

Quantum gases of atoms provide an ideal system for 
studying multi-component superfluids. At first sight, the 
rotational properties of a multi-component gas may look 
like a trivial generalization of the case of a single compo- 
nent. However, as long as the different components inter- 
act and exchange angular momentum, the extra degrees 
of freedom associated with the motion of each species is 
not at all a trivial effect. On the contrary, this coupled 
system may demonstrate some very different phenomena, 
see, e.g., Refs. @, H, [J]. Several experimental and theo- 
retical studies have been per formed on this problem, see, 
e.g., Refs. [j| [p, B , H, foL I lOL fill ], as well as the review 
article of Ref. [12| . 

In this Letter, the rotational properties of a super- 
fluid that consists of two distinguishable components are 
examined. Three new and surprising conclusions result 
from our study: 

Firstly, under appropriate conditions, one may achieve 
vortex states of multiple quantization. It is important 
to note that these states result from the interaction be- 
tween the different species, and not from the functional 
form of the external confinement. It is well known from 
older studies of single-component gases, that any external 
potential that increases more rapidly than quadratically 
gives rise to vortex states of multiple quantization, for 
sufficiently weak interactions [l3l |; on the contrary, in a 
harmonic potential, the vortex states are always singly- 
quantized. In the present study, vortex states of multiple 
quantization result purely because of the interaction be- 
tween the different components, even in a harmonic ex- 
ternal potential. Therefore, our study may serve as an 
alternative way to achieve such states [Tij ]. 

Secondly, our simulations indicate that multi- 
component gases do not support persistent currents, in 
agreement with older studies of homogeneous superflu- 
ids [3| . Essentially, the energy barrier that separates the 
(metastable) state with circulation/flow from the non- 
rotating state, is absent in this the numerical 



results, as well as the intuitive arguments presented be- 
low, suggest. 

Finally, we investigate the structure of the lowest state 
of the gas, in the range of the total angular momentum L 
between zero and iV m i n = min(N a, Nb), where Na and 
Nb are the populations of the two species labelled as A 
and B. In this range of L, only the single-particle states 
with m = and m = 1 are macroscopically occupied, as 
derived in Ref. fl5j | within the approximation of the low- 
est Landau level of weak interactions. Remarkably, our 
numerical simulations within the mean-field approxima- 
tion, which go well beyond the limit of weak interactions, 
show that this result is more general. 

For simplicity we assume equal masses for the atoms 
of the two components, Ma — Mb — M. Also, we 
model the elastic collisions between the atoms by a con- 
tact potential, with equal scattering lengths for colli- 
sions between the same species and different species, 
o-aa = clbb = clab = a (except in Fig. 4). Our results 
are not sensitive to the above equality and hold even if 
ciaa ~ a BB ~ o-ab, as in Rubidium, for example. For the 
atom populations we assume Na. ^ Nb, but Na/Nb ^ 1 
(without loss of generality). The trapping potential is 
assumed to be harmonic, V CK t(r) — M(u 2 p 2 + uj 2 z 2 )/2. 
Our Hamiltonian is thus 
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where Uq = 4ttTl a/M. We consider rotation around the 
z axis, and also assume that hu> z 3> hu>, and Tiuj z 3> no /To, 
where hq is the typical atom density. With these assump- 
tions, our problem becomes effectively two-dimensional, 
as the atoms reside in the lowest harmonic oscillator state 
along the axis of rotation. Thus, there are only two quan- 
tum numbers that characterize the motion of the atoms, 
the number of radial nodes n, and the quantum number 
m associated with the angular momentum. The corre- 
sponding eigenstates of the harmonic oscillator in two 
dimensions are labelled as <fr n ,m- 

Within the mean-field approximation, the energy of 
the gas in the rest frame is 
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where "J^ and & r e the order parameters of the two 
components. By considering variations in ty* A and ^f* B , 
we get the two coupled Gross-Pitaevskii-like equations, 
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where /iyi and /^s is the chemical potential of each com- 
ponent. We use the method of relaxation [l6| to minimize 
the energy of Eq. Q in the rotating frame, E' = E — Lf2, 
where fi is its angular velocity. 

For the diagonalization of the many-body Hamilto- 
nian, we further assume weak interactions, hqUq -C hu>, 
and work within the subspace of the states of the lowest 
Landau level, with n = 0. This condition is not neces- 
sary, however it allows us to consider a relatively larger 
number of atoms and higher values of the angular mo- 
mentum. We consider all the Fock states which are eigen- 
states of the number operators Na, Nb of each species, 
and of the operator of the total angular momentum L, 
and diagonalize the resulting matrix. 

Combination of the mean-field approximation and of 
numerical diagonalization of the many-body Hamiltonian 
allows us to examine both limits of weak as well as strong 
interactions. For obvious reasons we use the diagonal- 
ization in the limit of weak interactions, and the mean- 
field approximation (mostly) in the limit of strong in- 
teractions. The interaction energy is measured in units 
of w = U J\$ a Ax,y)\ 4 \Mz)\ 4 d 3 r = {2/7r) 1 / 2 hcja/a z , 
where (f>o(z) is the lowest state of the oscillator potential 
along the z axis, and a z = (h/AIujz) 1 ^ 2 is the oscilla- 
tor length along this axis. For convenience we introduce 
the dimensionless constant 7 = Nvo/Tiuj — y // 2/irNa/a z , 
with N = Na + Nb being the total number of atoms, 
which measures the strength of the interaction. 

We first study the limit of weak coupling, 7CI, and 
use numerical diagonalization. Considering Na = 4 and 
Nb = 16 atoms, we use the conditional probability dis- 
tributions to plot the density of the two components, 
for L = 4,16,28, and 32, as shown in Fig. 1. When 
L = 4 = N A , and L = 16 = Nb, the component whose 
population is equal to L forms a vortex state at the center 
of the trap, while the other component does not rotate, 
residing in the core of the vortex. This is a so-called 
"coreless" vortex state 

HHH2I- As L increases beyond 
L = Nb = 16, a second vortex enters component £?, and 
for L = 2Nb = 32, this merges with the other vortex to 
form a doubly-quantized vortex state. For this value of 
L = 32, the smaller component A does not carry any an- 
gular momentum (apart from corrections of order 1/N). 
The fact that this is indeed a doubly-quantized vortex 
state, is confirmed by the occupancy of the single-particle 
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FIG. 1: The conditional probability distribution of the two 
components, with Na = 4 (higher row), and Nb = 16 (lower 
row). Each graph extends between — 2.4<2q and 2.4ao in 
both directions. The reference point is located at (x, y) = 
(1.25ao,0) in the lower graphs (B component). The angular 
momentum L increases from left to right, L = 4(= Na), 16(= 
N b ), 28, and 32 (= 2N B ). 



states. By increasing Na, Nb, and L — 2Nb proportion- 
ally, we observe that the occupancy of the single-particle 
state with m = 2 of component B approaches unity, while 
the occupancy of all the other states are at most of order 
1/Ng. The same happens for the single-particle state 
with m = of the non-rotating component A. 

A similar situation emerges for the case of stronger cou- 
pling, 7 = 50, where we have minimized the mean-field 
energy of Eq. @ in the rotating frame (in the absence of 
rotation the two clouds do not phase separate). For ex- 
ample, we get convergent solutions, shown in Fig. 2, for 
Nb/Na = 2.777, and (i): L A = N A ,L B = 0, for CI/oj = 
0.35 (top left), (ii) La — 0,Lb = Nb, for Q/uj = 0.45 
(top middle), (iii) L A = 0.755N A ,L B = 1.171iV B , for 
n/w = 0.555 (top right), (iv) L A = 0,L B = 2N B , for 
fi/w = 0.60 (bottom left), (v) L A = 0.876N A ,L B = 
2.057 N B , for fi/w = 0.69 (bottom middle), and (vii) 
L A = 0,L B = 3N B , for Q/uj = 0.73 (bottom right). 
Here, La and Lb is the angular momentum of each com- 
ponent, with L = La + Lb- Again, when L — 2Nb, and 
L = 3Nb the phase plots show clearly a doubly-quantized 
and a triply-quantized vortex state in component B, and 
a non-rotating cloud in component A. 

The picture that appears from these calculations is 
intriguing: as f2 increases, a multiply-quantized vortex 
state of multiplicity k splits into n singly-quantized ones, 
and on the same time, one more singly-quantized vortex 
state enters the cloud from infinity. Eventually all these 
vortices merge into a multiply-quantized one of multiplic- 
ity equal to k + 1. Figure 2 shows the above results for 
various values of f2. 

The mechanical stability of states which involve the 
gradual entry of the vortices from the periphery of the 
cloud is novel. This behavior is absent in one-component 
systems, in both harmonic, as well as anharmonic trap- 
ping potentials. In one component gases, only vortex 
phases of given rotational symmetry are mechanically 
stable [H, In the present problem, the mechanical 
stability of states with no rotational symmetry (shown 
in Fig. 2) is a consequence of the non-negative curva- 
ture of the dispersion relation (i.e., of the total energy) 
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FIG. 2: The density (higher graphs of each panel) and the 
phase (lower graphs of each panel) of the order parameters 
(left graphs of each panel) and ^b (right graphs of each 
panel), with Nb/Na = 2.777 and a coupling 7 = 50. Each 
graph extends between — 4.41ao and 4.41ao in both directions. 
The values of the angular momentum per atom and of fl in 
each panel are given in the text. 
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£{L). This observation also connects with the (absence 
of) metastable, persistent currents (i.e., the second main 
result of our study), which we present below. 

In Ref. [l5[ we have given a simple argument for the 
presence of vortex states of multiple quantization within 
the mean-field approximation. At least when the ra- 
tio between Na and Nb is of order unity (but Na ^ 
Nb), there are self-consistent solutions of Eqs. ([3]) of vor- 
tex states of multiple quantization. Within these solu- 
tions, the smaller component (say component A), does 
not rotate, providing an "effective" external potential 
Kff,s(r) = K x t(r) + C/on^(r) for the other one (compo- 
nent B), which is anharmonic close to the center of the 
trap. This effectively anharmonic potential is responsible 
for the multiple quantization of the vortex states. There- 
fore, we conclude that for a relatively small population 
imbalance, the "coreless vortices" are vortices of multiple 
quantization. 

The second aspect of our study is the absence of 
metastable currents (in the laboratory frame, for fi = 0). 
A convenient and physically-transparent way to think 
about persistent currents is that they correspond to 
metastable minima in the dispersion relation £ = £(L) 
[20| . A non-negative curvature of £ (L) for all values of 
L implies the absence of metastability. For all the cou- 
plings we have examined, both within the numerical diag- 
onalization, as well as the mean-field approximation, we 
have found a non-negative second derivative of the dis- 
persion relation. Figure 3 shows La/Na, Lb/Nb, and 
L/N versus f2, for 7 — 50. These curves are calculated 
by minimizing the energy £{L) in the rotating frame for 
a fixed Q, and plotting the angular momentum per par- 
ticle of the corresponding state for the given rotational 
frequency. 

Again, our argument for the effective anharmonic po- 
tential is consistent with this positive curvature. Let us 



FIG. 3: Higher graph: The angular momentum La/Na 
(crosses) and Lb/Nb (dots), as function of Q. Lower graph: 
L/N as function of fi. All curves result from the minimiza- 
tion of the energy in the rotating frame, within the mean-field 
approximation, for 7 = 50. 



consider for simplicity claa = o-bb = 0, and clab 7^ 0. 
Then, the problem of solving Eqs. (3) becomes essen- 
tially a (coupled) eigenvalue problem. If Eo, m are the 
(lowest) eigenvalues of the effective (anharmonic) po- 
tential felt by the rotating component for a given an- 
gular momentum mTi, then d 2 EQ_ m /dm 2 is always posi- 
tive. For example, if one considers a weakly-anharmonic 
effective potential, V cS (p) = Mui 2 p 2 [l + \(p/a a ) 2k }/2, 
where k = 1, 2, . . . is a positive integer, ao = (fi/Mw) 1 / 2 
is the oscillator length, and < A < 1 is a small di- 
mensionless constant, according to perturbation theory, 
Eo,m — Tiw\m\ + A(|m| + 1) . . . (|m| + k)/2, which clearly 
has a positive curvature. 

One may gain some physical insight into the absence 
of persistent currents by understanding the difference be- 
tween a gas with one and two components. In the case 
of a single component, for sufficiently strong (and repul- 
sive) interactions, an energy barrier that separates the 
state with circulation from the vortex-free state may de- 
velop. In the simplest model where the atoms rotate in 
a toroidal trap, in order for them to get rid of the circu- 
lation, they have to form a node in their density, which 
costs interaction energy, and this creates the energy bar- 
rier [2(J, HH- On the other hand, in the presence of a 
second component, this node may be filled with atoms of 
the other species, and therefore the system may get rid 
of the circulation with no energy expense. This physical 
picture is also supported by the density plots in Figs. 2 
and 4. For example, in the case of coreless vortices, the 
core of the vortex is filled by the other (non-rotating) 
component fl2| . More generally, the density minima of 
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FIG. 4: The density (higher graphs) and the phase (lower 
graphs) of the order parameters (left graphs of each panel) 
and (right graphs of each panel), with Nb/Na = 2.777. 
Here Q/lj = 0.6 and 7 = 50. In the left panel La = 0, 
and Lb/Nb =2. In the right panel, the scattering length 
a,BB is twice as large as in the left panel, asB — 2a. In this 
case, La/Na = 0.05, and Lb/Nb = 1.936. All graphs extend 
between — 4.41ao and 4.41ao. 



the one component coincide, roughly speaking, with the 
density maxima of the other component, resulting in a 
total density n to t = \^a\ 2 + \^b\ 2 which does not have 
any local minima or nodes. 

Our third result is based on the mean-field approxima- 
tion. For < L < N min , where iV min = min(A^A, Nb), 
the only components of the order parameters ^>a and f ^ 
are the single-particle states with m = and m = 1, i.e., 

*A = J~] Cn,0$n,0 + Cn,i<i>n,i, 
n 

*B = ^d n) o$„,o+dn,l$n,l> (4) 

n 

where cv^o, Cn.,1, d ni o and d nj i are functions of L and of 
the coupling. The numerical simulations that we perform 
within a range of couplings 7 < 50 that extend well be- 
yond the lowest-Landau level approximation, reveal this 
very simple structure for the lowest state of both com- 



ponents. Also, the corresponding dispersion relation is 
numerically very close to a parabola, as in the case of 
weak interactions Again, one may attribute these 

facts to the effective potential that arises from the inter- 
action between the two species (l5j . 

In the studies that have examined a single-component 
gas in an external anharmonic potential, it has been 
shown that as the strength of the interaction increases, 
there is a phase transition from the phase of multiple 
quantization to the phase of single quantization [19j | . In 
the present case the situation is more complex, since the 
effective anharmonic potential is generated by the in- 
teraction between the two species as a result of a self- 
consistent solution. Still, a similar phase transition takes 
place here, as, for example, one keeps the scattering 
lengths aAA and ciab fixed, and increases clbb that cor- 
responds to the rotating component. Figure 4 shows the 
density and the phase of both species, for aAA = clbb = 
o-ab = a (left panel), and clbb = 2a a a = 2a ab = 2a 
(right panel). Component B undergoes a phase transi- 
tion from a doubly-quantized vortex state to two singly- 
quantized vortices. 

To conclude, mixtures of bosons demonstrate numer- 
ous novel superfluid properties and provide a model sys- 
tem for studying them. Here we have given a flavor of 
the richness of this problem. Many of the results pre- 
sented in our study are worth investigating further, as, 
for example, one changes the ratio of the populations, the 
coupling constant between the same and different species, 
or the masses. 
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